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Helical ribbons arise in many biological and engineered systems, often driven by 
anisotropic surface stress, residual strain, and geometric or elastic mismatch be- 
tween layers of a laminated composite. A full mathematical analysis is developed 
to analytically predict the equilibrium deformed helical shape of an initially flat, 
straight ribbon, with prescribed magnitudes and orientations of the principal cur- 
vatures when subjected to arbitrary surface stress and/or internal residual strain 
distribution. The helix angle, radius, axis and chirality of the deformed helical 
ribbons are predicted with a comprehensive, three-dimensional model that incor- 
porates elasticity, differential geometry, and variational principles. In general, the 
mechanical anisotropy (e.g., in surface/external stress, residual strain or elastic 
modulus) will lead to spontaneous, three-dimensional helical deformations. Ring 
shapes are formed when the principle axes of deformation coincide with the geo- 
metric axes of the ribbon. The transition from cylindrical helical ribbons to purely 
twisted ribbons, or tubular structures is controlled by tuning relevant geometric pa- 
rameters. Analytic, closed-form predictions of the ribbon shapes are validated with 
simple, table-top experiments. This theoretical approach represent a tool to inform 
the design of materials and systems in order to achieve desired helical geometries 
on demand. 

Keywords: Helices, Elasticity theory, Surface stress, Nanohelices, Tubules 

1. Introduction 

Helical structures arise in a variety of natural systems, such as plant tendrils 
(Chouaieb et al. 2006; Gcrbode et al. 2012) , chiral seed pods (Armon et al. 2011; 

f To whom correspondence should be addressed. E-mail: chen.z@seas.wustl.edu. 
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Forterre & Dumais 2011) , stork's bill awn (Abraham et al. 2011) and twisted guts 
(Savin et al. 2011; Wyczalkowski et al. 2012) , as well as nanoengineered polymer 
lamellae (Wang et al. 2011) , twist-nematic-elastomer films (Sawa et al. 2011) , 
ZnO nanohclices (Kong et al. 2003; Majidi et al. 2010) , graphene nanoribbons 
and nanotubes (Kit et al. 2012; Cranford and Buehler 2011) . They have since at- 
tracted scientific attention for their potential role in nanoelectromechanical systems 
(NEMS) (Bunch et al. 2007), drug delivery and biological/chemical sensing(Cui et 
al. 2001; Smith et al. 2001;Zastavker et al. 1999), magnetic field detection (Huang 
et al. 2005), optoelectronics (Hwang et al. 2008), and microrobotics (Abbott et 
al. 2009). With recent advancements in nanotechnology, physicists and engineers 
can now grow helical nanoribbons (Kong and Wang 2003) through a "bottom-up" 
approach and have also begun exploring ways to fabricate helical nanoribbons of 
controllable morphology (Bell et al. 2007; Prinz et al. 2000; Zhang et al. 2005) in a 
"top-down" manner (Cho et al. 2006). Motivated by potential applications in pro- 
teomics and biological surface engineering, scientists have also explored methods to 
control the morphology of helical nanoribbons through hierarchical self-assembly 
of amphiphilic molecules (Zhang et al. 2002). Helical nanoribbons may also play a 
key role in emerging technologies like programmable matter, which adopt a broad 
range of geometries and functionalities in a controlled and reversible manner. 

In this manuscript, we examine the mechanics and morphology of thin, effec- 
tively two-dimensional elastic helical ribbons embedded in within a theoretical 
framework that accounts for elastic bending along both the length and width direc- 
tions (Chen et al. 2011) . Previous efforts have largely focused on one-dimensional 
representations for use in molecular dynamics simulations (Lee et al. 2011; Vnai 
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and Zakrzewska 2004) and continuum calculations (Chouaieb et al. 2006; Healey 
and Mehta 2005, Helfrich 1986; Fuhrhop and Helfrich 1993; Kirchhoff 1859) that 
treat the helix as a space curve. While this approach may be appropriate for sys- 
tems like DNA (Bryant et al. 2003; Snir and Kamien 2005) , a-helices (Shepherd 
et al. 2009) , and carbon nanocoils (Bandaru et al. 2007) , there exist an important 
class of unique helical morphologies in nature and engineering (Chung et al. 1993; 
Giomi and Mahadevan 2010; Fuhrhop and Helfrich 1993; Lugomer and Fukumoto 
2010; Selinger et al. 2004; Srivastava et al. 2010) that can only be captured with a 
two-dimensional representation that treats the helix as a ribbon of finite width. The 
helical ribbon model, presented here, represents a two-dimensional generalization 
of the classical Stoney/Timoshenko model used to study planar bending (Stoney 
1909; Suo et al. 1999; Timoshenko 1925; Zang et al. 2007). 

In many micro- and nano-scale systems, the mechanics of helices is governed by 
surface stress and other residual stresses that are comparable to the elastic restoring 
forces of bulk deformation (Gurtin and Murdoch 1975; He and Lilley 2008; Li et 
al. 2010; Lotz and Cheng 2005; Miller and Shcnoy 2000; Lachut and Sader 2007) . 
These stresses may arise from surface defects or imbalances in atomic coordination 
number, adsorbates, and/or surface reconstruction (Yi and Duan 2009). The shape 
of the helix may be controlled not only by the magnitude of the surface stress but 
also the anisotropy of both the residual stress field and ribbon elasticity (Chen et 
al. 2011; Lotz and Cheng 2005; Wang et al. 2008; Ye et al. 2009; Zhang et al. 2005; 
Zhang et al. 2006; Wang et al. 2011; Wang et al. 2012) For example, helicity may 
be introduced by the orientation of the in-plane principle components of either the 
elastic stiffness or surface stress, or a combination of both. The equilibrium shape 
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may, in turn, be influenced by interactions with gas phase composition (Bandaru 
et al. 2007) , surface chemistry (Randhawa et al. 2010) , temperature (Timoshcnko 
1925) , solution pH (Zhang et al. 2010) , or other environmental factors that might 
affect surface stresses. For ribbon laminates composed of two or more layers, helicity 
can also be controlled by lattice mismatch or pre-stretching individual layers (Zhang 
et al. 2005) . Although the physical origins of these stress asymmetries are different, 
we will demonstrate that the mechanics of each system represents special cases of 
a more general mathematical description. 

Typically, ribbon helicity arises because at least one of the mechanical elements 
(e.g., surface stress (Wang et al. 2008) , external stress (Pokroy et al. 2009) , residual 
strain, or elastic modulus (Zhang et al. 2005) is anisotropic and the principle axes 
of curvature do not coincide with the principle geometric axes (length and width) 
of the ribbon. For quaternary sterol solutions such as model bile, anisotropic stress 
may be caused by a mismatch in molecular packing between constituent layers 
(Zastavker et al. 1999) , while for nano-engineered helices, the deformation may 
be driven by epitaxial strains (Hwang et al. 2008) . Note that helical deforma- 
tions represent large deflections of thin plates and is therefore a complex, nonlinear 
geometric problem that is usually difficult to solve even in very simple cases (Lan- 
dau and Lifshitz 1986) . However, there has been increasing interest in studying 
morphologically distinct deformations such as purely twisted ribbons (with saddle 
curvature) (Liu et al. 2009; Oda et al. 1999) , cylindrical helical shapes (Zhang et 
al. 2005) and ring configurations (Cho et al. 2010; Mei et al. 2009) , as well as in 
the connections between them (Bellesia et al. 2008; Ghafouri and Bruinsma 2005; 
Iglic et al. 2005; Selinger et al. 2001; Selinger et al. 2004) . Although tradition- 
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ally each geometrically distinct object (Iglic et al. 2005; Selinger et al. 2001) has 
been studied separately, recent studies (Bellesia et al. 2008; Ghafouri and Bruinsma 
2005; Selinger et al. 2004) have explored the transitions between purely twisted rib- 
bons and cylindrical helical shapes. In the present work, we show that all of the 
above-mentioned shapes, in fact, belong to the same geometric class, governed by 
only three tunable geometric parameters. We demonstrate that transitions between 
these different shapes occurs smoothly as the geometric parameters are varied. 

In this paper, we first employ concepts from differential geometry to establish 
a mathematical description of a ribbon deformed with a given set of principle cur- 
vatures (magnitudes and orientations) , and show that the centerlinc deforms into a 
circular helix with constant radius and pitch. The as-deformed ribbon morphology 
represents a class of two-dimensional surfaces controlled by three independent geo- 
metric parameters, which, when properly tuned, give rise to different shapes such as 
cylindrical helical ribbons, purely twisted ribbons and ring shapes. Then, we employ 
continuum elasticity theory and stationarity principles to quantitatively link princi- 
ple curvatures and surface stresses or residual strains, providing a capability to pre- 
dict the morphology of stress/strain- induced helical ribbons. Closed-form analytic 
predictions are obtained and validated with simple, table-top experiments where 
layers of pre-stretched elastic sheets are bonded together to form laminated rib- 
bons. Good agreement is achieved between our theoretical analysis, which uniquely 
predicts shapes that cannot be explained with classical rod or plate theories, and 
experiments with no adjustable parameters. Also, because the ribbon bends in 
two directions with non-zero Gauss curvature, our analysis is beyond the scope of 
the classical Stoney formulation of planar bending of ribbons under surface stress 
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Figure 1. Illustration of a helical ribbon. The vectors d x and d y are oriented along the 
length- and width-wise axes of the ribbon, respectively. The bases ri and r-2 correspond 
to the principle directions of curvature. 

(Stoney 1909) . This paper provides detailed derivations of the geometrical and 
continuum elasticity model along with explicit results and discussions, and hence 
greatly expands on the work in a very brief companion manuscript on this topic 
published previously by the authors (Chen et al. 2011) . 

2. Geometry of a helical ribbon 

A ribbon is defined as an elastic strip of length L, width w <C L, and thickness 
H -C w (Chen et al. 2011) . The cross-section of the ribbon is rectangular and the 
principle geometric axes of the ribbon are oriented along the ribbon length (d^), 
width (d y ), and thickness (d z ) directions. These directions form an orthonormal 
triad of vectors {d x ,d y ,d z } that rotate with the ribbon as it coils in a three- 
dimensional Cartesian frame {E x , E y , E z }. This representation originates from the 
directed rod theory in Green & Laws (2006) and ORcilly (1998) 

In planar bending (Stoney 1909) , curvature is restricted to the (length-wise) 
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direction. In the general case of helical deformation, illustrated in Fig. [TJ curvature 
will occur along principle directions 

ri = d x cos <j) — d y sin0 and r 2 — d x sin <fi + d y cos cf> , (2-1) 

where <f> represents the orientation of the principle curvatures with respect to the 
principle geometric axes and K\ and k 2 denote the principle curvatures along the 
directors ri and r 2 , respectively. Points along the centerline of the ribbon follow a 
space curve 

P = P (s) = X Q (s)-E x + Y (s)-E y + Z (s)-E z , (2.2) 
where s G {0, L} is the convecting arclength. In the following, we derive a rela- 
tionship between the ribbon curvatures and orientation {ki,K2,4>} and the pitch, 
radius, axis, and chirality of the helical space curve P = Po(s). In doing so, it is 
convenient to define the unit normal N = d z = x d y — ri x r2. By definition, 

— — - = — Nki cos 4> an( l —j — = — NK2 sin <f> . (2-3) 
as as 

Lastly, the identities dP/ds = d x and N = ri x r 2 together imply 

dP . , , , dN , . , ,„ 

— — = sm 0ri + cos 0r2 and — — = m cos (pri + K2 sm cpr2 ■ (2-4) 
as as 

(a) Derivation of the helix centerline coordinates 

Next, we derive a differential equation for N(s) and solve it explicitly. To this 
end, from Eqs. ([2~3]l and (|2~4")l it directly follows that 
d 2 N 

— = - (k? cos 2 <j> + K 2 2 sin 2 <P) N . (2.5) 

This second-order differential equation is subject to the initial conditions N(0) = E z 
and 

= Kl cos0ri(O) + K 2 sin0r 2 (O), (2.6) 

as 
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n (0) = cos <fi~E x - sin <p~E y and r 2 (0) = sin <j>E x + cos <j>Ey . (2.7) 
Solving for N(s) yields 



N = ^— sin as^j E x — ^— sin as^ E y + (cos as) E 2 



where 



(2.8) 



a = y k 2 cos 2 (f> + k 2 , sin 2 0, (2-9) 
/3 = ki cos 2 4> + k 2 sin 2 0, (2.10) 

r = («i — k 2 ) sin</>cos(/) . (2-11) 



and 



In passing, we note that a 2 = fP + r 2 . Next, integrating Eqs. (|2.3[) for the initial 
conditions expressed in Eq. (|2.7I) yields 

ri = icos</> + -^Kifcosas — 1) cos^l E x 
l" 2 J 

+ |— sin0 -ki(cos as — 1) cos j E y 

- | ^- sin as cos </)})E 2 (2.12) 

r 2 = { sin 6+ —r K2 (cos as — 1) sin 0> E x 
I " J 

+ |cos 4> — ^-K2(cos as — 1) sin</>j E y 
-|^sinassin0|E z . (2.13) 

Lastly, integrating the expression for dP / ds in Eq. (|2.4I) implies 



Xq(s) — s t (as — sin as) 

a J 

3^b (s) = — r- (as — sin as) 
a J 

Z (s) = A( cosas _i). (2.14) 
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As a consistency check, we see that K\ — k 2 — 0, Xq(s) = s while Yq(s) — Zq(s) = 
describes a straight ribbon oriented along the E x direction. 

(b) Geometric properties of the helical space curve P 

The space curve P derived in Sec. m represents a helix with curvature j3, torsion 
r, axis M, helix angle $ (the angle between a tangent vector of the helix and the 
plane normal to the helix axis), radius R, and pitch D. The corresponding Frenet- 
Serret frame {T,N, B} (Tu and Ou-Yang 2008) satisfies the following differential 
equations: 

^=- J 8N(a) ^=/?T(s)+rB(s) ^=-tN(«). (2.15) 

Here, T = is the tangent vector, 

T(s) = (l- ^-(1 - cos as) 1e x + ^(1 -cosas)E v - ( -sin as) E z , (2.16) 
I a z J a z I a J 

and B = T x N (Kamien 2002) is a bi-normal vector: 

— cosas)E x + < — ^(1 — cos as) + cos as 1 E„ + {— sin as} E z , (2.17) 
a 2 ya 2 J L a J 



B 



Now, the axis M is parallel to the vector formed by one complete repeating unit: 
P(27r/a) - P(0) = 2fJ {^E, + |E S }. Normalization yields M = ^E x + &E V , 
implying that, for a helix, T • M = r/a is a fixed value, which corresponds to the 
helix angle 

. fT\ ( (ki — K2) sin0cos</> 1 

$ = arcsin - = arctan <^ ^ ^ " — ^- \ . (2.18) 

v a) \n\ cos 2 </> + k 2 sin J 

Helical chirality is determined by the sign of the helix angle (or equivalently the 

torsion of the ribbon centerline), i.e., sgn($). Right handed helices correspond to 

$ > whereas left-handed helices correspond to $ < 0. Equation (|2.18p implies 
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that the chirality is determined by the signs of both («i — k%) and sin cos 0. As 
expected, when interchanging the values of K\ and K2, the helix shape transforms 
to its mirror image about the plane y = z with reversed handedness. Lastly, the 
helix has a radius 

or k\ cos 2 (j) + sin (j> 



and pitch (jjj) 



?7r-r 

D = ||P(27r/a)-P(0)|| = - r . (2.20) 



(c) Helical space curve: Special cases 

We now consider a few special helical shapes before describing the more general 
case. The ribbon forms a ring when either k,\ — k,%, (f> — 0, or (f> = w/2, which leads 
to helix angle of zero according to Eq. (|2.18|) . For example, when n\ = K2 [as shown 
in Fig. (He)], the helix angle is $ = and the helix reduces to a ring. A ring can also 
be obtained when <j> — or 7r/2. Moreover, when the Gauss curvature K — k\k,i is 
zero, a helical cylindrical ribbon is formed. For example, setting K2 — [cf. Fig.JSJa)] 
yields a right-handed helical ribbon that wraps around a cylinder of radius l/«i, 
while setting k\ — [cf. Fig. EJd)] leads to a left-handed helical ribbon wrapping 
around a cylinder of radius 1/k%. Another configuration of particular interest is 
the purely twisted ribbon, which forms when n\ cos 2 <fi + K2 sin 2 = implying 
that the helix angle is $ = ir/2 and radius R = 0. A specific example is shown in 
Fig-IH/) for the case K2 = — K\ (i.e., the mean curvature H = (kj + K2)/2 = 0) and 

f An alternative method for deriving R and D is to examine Eq. 112. 1511 directly. Since the 
Prenet-Serret triad associated with the centerline satisfies Eq. A2.15I I, with constant curvature 
P = K\ cos 2 <f> + «2 sin 2 ij> an< l torsion r = (fti — «2 ) sin <j> cos <j>, the centerline is indeed a circular 
helix with constant radius R = /3/(/3 2 + t 2 ) = /3/a 2 and pitch _D = 27rr/a 2 . 
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Figure 2. Spectrum of representative ribbon morphologies as relative values of Ki and k,2 
are tuned (<f> = 7r/4 in all these cases). As described in Section[2] the ribbon originally lies 
along E x direction in global coordinate system (which forms an angle <f) from r?), in the 
plane spanned by rj and r° where r° = cos^E^ — sin^E^ and r% = sin^E^ + cos(j!>E y . 
The values of principal curvatures Ki and K2 are as follows: (a) ki = 1, «a = 0; (b) 
Ki = 1,K2 = 1/2 ; (c) «i = K2 = 1; (d) kj = 0, K2 = 1; (e) «i = 1/2, = 1 ; (f) 
«i = — 1, «2 = 1 . 

= 7r/4. Here, the ribbon undergoes a pure twist deformation (i.e., the centerline 
remains straight). In other words, by tuning the values of Ki, «2 (or equivalently the 
mean and Gauss curvatures) and <$> (in the above cases, <f> is fixed to be 7r/4, but in 
general, 4> is another tunable geometric parameter) , enables one to generate helical 
shapes of any kind. More specifically, transitions between rings, cylindrical helical 
ribbons, and purely twisted ribbons can be achieved by tuning the set of three 
independent geometric parameters, namely the two principal curvatures {ki,k^) 
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and the misalignment <p between principal curvature and geometric axes. For given 
«i, K2, and <p, the ribbon morphologies were generated by the procedure outlined 
in the Appendix. 

3. Continuum elasticity theory and ribbon morphology 

Having established a quantitative link between prescribed ribbon principal curva- 
tures and ribbon morphology in Sec. [H we now focus on the role of surface stress 
in controlling the principal curvatures (and hence morphology). To this end, we de- 
termine the equilibrium shape of the helical ribbon by minimizing its elastic strain 
energy, 

n = f- : 7U=-i?/2 + f + : 7U=ff/2 + / ' \l ■ C : 7 dz . (3.1) 

J-H/2 1 

Here, 7 is the strain tensor for linear elastic deformation, C is the fourth-order 
elastic stiffness tensor, and {f + ,f - } are the surface stresses. The coordinate z £ 
{—H/2, H/2} is defined through the thickness of the ribbon and f + and f _ corre- 
sponding to the surfaces at z = H/2 and z = —H/2, respectively. 

(a) Kinematics and stationarity principle 

As in classical plate theory, the strain tensor 7 is obtained by superimposing 
strains induced by elastic bending (7^) and an in-plane uniform "membrane" strain 
(7m)- We also include strain along d z in order to allow for plane stress compatibility 
(7 Z ). In general, the ribbon may also be subject to internal residual strains that 
do not correspond to surface stress (7 ). Combining these contributions leads to a 
strain tensor of the form 

7 = 76 + 7m + 7 Z + 7o , ( 3 - 2 ) 
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where 

7 b = zkiTi ® ri + zK 2 r 2 ® r 2 

7 m = £xxd x (8>d a: +£ X j / (d x ®d ?/ + dj / (8)d x )+£ TO dy (8)^ 

7* = (£zz + qz)d z <g> d z 

7o = 7 ° x d x ® d x + 7 ° y (d x ® d y + d y ® d x ) + 7° y d y ® d y + 7 2 ° z d z ® d^.3) 

Here, g denotes the gradient of the strain component along d z required for plane 
stress compatibility. 

The strain tensor may also be expressed as 7 = 7yd j ® d^ (i, j e {x, y,z}) with 
components 

7xx = £xx + z(m cos 2 + k 2 sin 2 <f>) + -f xx {z) 

Ixy = lyx = £xy + z(k 2 - Kl ) sin COS + 7° y (z) 

Ivv = £ yy + z («i sin2 ^ + «2 cos 2 <j>) + l° yy {z) 

Izz = e zz + qz + ~t° zz {z). (3.4) 

Let <j) + and 0~ denote the orientation of the principle axes of the surface stress 
on the top (z = if/2) and bottom (z = — if/2) surfaces, respectively. Unlike cj>, 
which is an unknown geometric parameter, both cf) + and cf)~ represent fixed values, 
prescribed by nature or in manufacturing. For example, they could coincide with the 
principle axes of the surface stress tensor, f. On the top and bottom (z = ±if/2) 
surfaces, f adopts the following forms: f ± = f^ef ® ef + f2 e 2 ® e t> where 

= cos0 ± d 2: — sin(/) ± d 2/ and = sin^ ± d 2; + cos^ ± d !/ . 

At equilibrium, the variational principle dictates that II must be stationary with 
respect to the variations of all the unknown parameters K\, k 2 , q, e XI , e yy , e xy , e zz 
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and 4>. In other words, these values correspond to the solutions to the following 
set of linearly independent equations: dH/dx — 0, where \ represents any one of 
the eight unknowns. It is worth noting that although here we choose to write the 
strain tensor in the {d a; ,dj ) ,d z } coordinate system, it could also be expressed in 
other orthonormal coordinate systems, such as {ri, T2, d z } (which we will adopt in 
the example that follows). The total potential energy density, however, remains the 
same regardless of which orthonormal coordinate system one chooses. 



In general, the system of eight equations dH/d\ = for the eight unknowns 
{x} can only be obtained numerically. However, closed-form algebraic solutions are 
possible when the ribbon is treated as homogenous, elastically isotropic, and subject 
to surface stress on only one surface. Without loss of generality, we consider surface 
stress acting on the bottom surface {z = —H/2) only: f~ = /•j - ej"<g>ej" + - /2~e^"<g>ej, 
and the angle from d^ to (clockwise) is <p~ . As illustrated in Fig. [TJ the strain 
tensor is diagonalized in the orthonormal coordinate system {ri,r2,d z }, so that 
the only nonzero strain components are 711 = en + zk\, 722 = £22 + zk 2 and 
733 = £33 + Q z - Thus, the total potential energy per unit area is 



(b) Analytical solutions 



n 




+ 



EH { v{e 11 +e 2 2 + £33? (4i + 4> + £33) 
2 \ (1 - 2i/)(l + v) l + i/ 



} 



(3.5) 
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where E and v denote the Young's modulus and Poisson's ratio of the ribbon, 
respectively. 

Applying the variational principles implies that the principle directions of the 
curvature and surface stress coincide, i.e., <p — 4>~ ■, an d that 

6(/r - v£) K _ 6(/ 2 --^/r) 
Kl ~ Em K2 ~ 

6K/r + / 2 _ ) . fr-vfc 



1 = FTFT^ £ H = - 



EH 2 11 EH 

e22 = _/2--^/f e33 = ^ML+in. (3.6) 
22 EH 33 EH V ' 

It is interesting to note that in the absence of Poisson coupling, the first two equa- 
tions are reduced to the classical Stoney formula (Stoncy 1909) in the two principle 
curvature directions: K\ = 6/f /EH 2 and k 2 = 6/^~ /EH 2 . 

For the more general case, where surface stresses occur on both surfaces, we 
can decouple the stretching and bending modes. First, we define the stretching 
component of the surface stress f s = f + + f and the bending component f * = 
f — f+, where f* represents the effective surface stress acting on the bottom surface 
(z = —H/2). Next, f s and f* can be diagonalized with respect to the principle axes 
(ui, u 2 ) and (n, r 2 ): f s = /fui <g> u x + ftu 2 <g> u 2 and f* = ftn ® r 1 + / 2 *r 2 <g> r 2 . 
Thus, by de-coupling the stretching and bending modes, and applying variational 
principles, the total strain tensor 7 = 7 S + 7/, can be obtained. 7 S is the strain 
component due to stretching 

ft - yfi _ ft - "ft » H.ft + .fi) ^ _ . 

ls = EH~ Ul 1 EH~ U2 U2 £ff z z ^ ^ 

and 7 b is the strain component caused by an effective surface stress on the bottom 
surface 

7 6 = Kizri <g) ri + K 2 zr 2 ® r 2 + qzd z ® d z . (3.8) 
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Figure 3. Fabrication of a helical ribbon. A strip of elastic, pressure sensitive adhesive is 
bonded to a pre-stretched sheet of latex rubber, with a misorientation angle <f>~ between 
the principal curvature and geometric axes. 

Thereby, the problem is reduced to an equivalent system with surface stress only 
on the bottom surface, with n x = 6(ft - vft)/EH 2 , n 2 = 6(/ 2 * - vf?)/EH 2 , and 
q = -&v{r i +t 2 )/EH\ 

4. Experimental validation and discussion 

We validated our theoretical predictions with a series of simple, table-top experi- 
ments in which a thin sheet of latex rubber was pre-stretched and bonded to an 
elastic strip of thicker, pressure-sensitive adhesive (Fig. [3]), similar to the physical 
models employed by the current and other researchers (Chen et al. 2011; Armon 
et al. 2011; Chen et al. 2012; Huang et al., 2012; Gerbode et al. 2012) . The latex 
sheets were purchased from Small Parts Inc. and the elastic strips were Acrylic, 
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Scotch Wall-Mounting Tape, produced by 3M. Their Young's moduli were 1.4MPa 
and 10.3MPa respectively, measured by an Instron™ 5848 Micro- Tester, and the 
thicknesses were 0.048cm and 0.1cm respectively. The Poisson's ratio of Acrylic is 
0.37 (Powers and Caddell 1972) and that of the latex sheets is 0.49 (Kaazempur- 
Mofrad et al. 2003) . Helical ribbons with different helix angle, radius, axis, chirality, 
and Gauss curvature were obtained by altering the magnitudes of the two princi- 
ple pre-stretches and their angle with respect to the centerline orientation of the 
bonded strip. These results were then compared with the theoretical predictions, 
as described in the captions to Figs. 0HH] In Figs. @J the predicted principal radii of 
curvatures are 19.6cm and 123cm respectively, consistent with the experimentally 
measured values 17.6cm and 112cm. In Fig. [5j the predicted principal radius of 
curvature is 27.0cm, in agreement with the experimental value 23.0cm. While for 
the purely twisted ribbon as shown in Fig. [51 the theoretical calculation and the 
experimentally measured value of the principal radii of curvature are ± 12.4cm and 
±13. 6cm respectively, again showing reasonable agreement. 

When the pre-stretches are anisotropic, helical ribbons can exhibit a broad range 
of helix angles, radii, and centerline orientations. As shown in Fig.|H the theoretical 
predictions for helix angle and ribbon shapes are in good agreement with the ex- 
perimental observations with no adjustable parameters. We note that experimental 
validation of the theory has also been demonstrated for the more general case in 
which both sides of the elastic strip arc bonded to differently oriented pre-stretched 
layers of latex rubber (Chen et al. 2011) 

It is important to emphasize that helicity in these systems arises from me- 
chanical anisotropy. Such anisotropy may originate from elastic anisotropy due to 
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Figure 4. A square piece of latex rubber (solid lines) is pre-stretched twice as much in 
the vertical direction than in the horizontal direction (to the dash-dot edges). An elastic 
adhesive sheet with no pre-strains is then bonded to the strained latex sheet and subse- 
quently cut into strips with long axes varying between 0° and 90° (at an interval of 15°). 
The relaxed ribbons are shown at the appropriate angles along with the corresponding 
theoretical predictions. 

molecular tilt or chiral interactions (Selinger et al. 2004) , or mechanical anisotropy 
in surface/external stresses, residual strains and elastic properties (Zhang et al. 
2005; Wang et al. 2008; Dai and Shen 2009; Pokroy et al. 2009) . In molecular dy- 
namics simulations, anisotropy can also arise from the asymmetry in the direction 
along which periodic boundary condition is applied (Lee et al. 2011) . The presence 
of at least one of these sources of anisotropy is required in order to break the ge- 
ometric symmetry and define the chirality (handedness) of the helical ribbon. For 
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single-layered ribbons, helicity can be introduced by the anisotropy of the in-plane 
principle components of either the elastic stiffness or surface stress tensors, or a 
combination of both. 

Although our current model only accounts for the effects of surface stress anisotropy, 
the same principle can be applied to treat other types of anisotropy that lead to 
helicity, such as helicity in ribbon laminates composed of two or multiple layers due 
to lattice mismatch or pre-strained individual layers (Zhang et al. 2005) . Alter- 
natively, if the pre-stretches are isotropic (biaxial tension), ring shapes will result, 
regardless of the different centerline orientations, as often observed in the literature 
(Mei et al. 2009; Cho et al. 2010) . Figure [5] shows good agreement between the ex- 
perimental results and theoretical predictions with such equal-biaxial pre-stretches 
on the bottom elastic sheet. 

Purely twisted ribbons represent yet another class of interesting morphologies 
(Oda et al. 1999; Selinger et al. 2001; Liu et al. 2009) . Such a deformation is outside 
the scope of developable ribbon deformations (Giomi and Mahadevan 2010) , since 
the deformation of a purely twisted ribbon is non-isometric (i.e., local distances 
between points are not always preserved). Since our theory takes into account all 
possible deformation modes (i.e., stretching, bending and twisting) of the ribbon, 
it has the capability of capturing pure twist. As predicted by the current theory, 
when Ki cos 2 <fi + Ki sin 2 = 0,$ = vr/2, a purely twisted ribbon will form. This 
was also verified in our experiments, as shown in Fig. [5] The rectangular bottom 
sheet (ABCiDi) and the top sheet (AB2C2D) are pre-stretched in the vertical and 
horizontal direction respectively so that they both become a square sheet (ABCD), 
before bonded with an unstrained elastic strip with <fi = 45°. According to our the- 
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Figure 5. A square piece of latex rubber (solid lines) is stretched equally in the vertical 
and horizontal direction (dash-dot edges). An unstrained elastic adhesive sheet is then 
bonded to the strained latex sheet and subsequently cut into ribbons with long axis varying 
between 0° and 90° (15° intervals). The released samples are shown at the appropriate 
angles along with the corresponding theoretical predictions. All of the ribbons adopt ring 
shapes with identical radii regardless of the orientation. 

ory, this yields K2 = — Ki (the same as shown in Fig. [5]) resulting in a purely twisted 
ribbon; the predicted morphology is in excellent agreement with the corresponding 
experimental result. 

Structural transitions among cylindrical helical ribbons, general helical ribbons, 
purely twisted ribbons, and cylinders/tubules are all captured within the current 
theoretical framework, as illustrated in Figs. [2] and [7] Similar transitions have been 
observed in bio-chemical systems such as cholesterol (Chung et al. 1993) , surfac- 
tants (Oda et al. 1999) and peptides (Bellesia et al. 2008; Pashuck and Stupp 2010) 
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Figure 6. A rectangular bottom sheet (ABCiDi) and a top sheet (AB2C2D) are 
pre-stretched in the vertical and horizontal direction, respectively, so that they both 
become a square sheet (ABCD), before bonding with an unstrained elastic strip with 
<j> = 45°. As predicted, the strip deforms into a purely twisted ribbon. 

, as well as in inorganic materials such as the transition metal dichalcogenide WS2 
(Iglic et al. 2005) . Understanding these transitions is key to programmable fab- 
rication of nano- and bio-devices, and possibly, a better understanding of certain 
diseases where morphological changes of helical ribbons and tubules are key (Chung 
et al. 1993) . More specifically, our kinematic model shows that when the Gauss 
curvature K is negative, the ribbon takes on a local saddle shape (the first two 
panels in Fig. [7]), and in particular, if n\ cos 2 <p + K2 sin 2 = 0, the ribbon becomes 
a purely twisted strip (the first panel in Fig. [7]) as investigated by (Selinger et al. 
2001) , (Ghafouri and Bruinsma 2005) , and (Iglic et al. 2005) 
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When the Gauss curvature vanishes, a helical ribbon with cylindrical curvature 
(like a ribbon bounding exactly about a cylinder so that the edges are straight) 
adopts the third configuration in Fig. [Jj However, for non-zero mean curvature and 
non-zero Gauss curvature, the general configuration is still helical, as shown in the 
second and fourth configurations of Fig. [7J For these ribbons, the edges are not 
exactly straight and could be either concave or convex depending on the sign of the 
Gauss curvature. This represents a class of general helical configurations that have 
been largely overlooked in the literature but may contribute to the understanding 
of the gradual conformational changes from twisted to helical ribbons that occur in 
a variety of physical systems (Fuhrhop and Helfrich 1993, Oda et al. 1999, Selinger 
ct al. 2004) . On the other hand, the last two configurations in Fig. [7J show the 
formation of "tubules" when the edges of the ribbon touch; such configurations can 
be obtained by either increasing ribbon width w at fixed misorientation <j) or by 
varying </> at fixed w. The as-generated tubules are similar to those formed in chiral 
liquid crystals due to molecular tilt modulation as illustrated in Fig. 6 of (Selinger 
et al. 2001) , although the origins of the mechanical anisotropics are presumably 
different. 

Historically, there have been limited theoretical attempts to explain the change 
between purely twisted ribbons with saddle curvature (K < 0) and helical cylin- 
drical geometries (K = 0) (Bellesia et al. 2008; Ghafouri and Bruinsma 2005; Iglic 
et al. 2005; Selinger et al. 2001; Selinger et al. 2004) . The existing models often 
attribute the morphological changes to the asymmetry in the chemical nature and 
the affinity of the two molecular surfaces to the solvent, or the force differences 
on two surfaces that lead to preferential bending or twisting. These asymmet- 
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Figure 7. Structural transition from twisted to helical ribbons and tubules with varying 
principal curvatures K2), misorientation cj>, and width w. 

ric conditions, when applied in molecular dynamics or Monte Carlo simulations, 
yield configurations consistent with experiments (Bellesia et al. 2008; Selinger et 
al. 2004) . Nevertheless, these approaches are computationally expensive, and find- 
ing an energy minimum often requires a starting configuration near the minimum 
energy state for fast convergence. In contrast, the current theory, by conveniently 
parameterizing those asymmetries through surface stress or elasticity, provides a 
continuum approach towards quantitatively understanding the associated struc- 
tural transitions induced by tuning the relevant geometric parameters. Of course, 
proper coarse-graining models need to be developed to establish the link between 
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the asymmetry /anisotropy in molecular interactions and the effective, anisotropic 
surface stresses that result in helicity. 

5. Conclusions 

A mathematical analysis has been developed to predict the formation of families 
of helices, rings, and twisting ribbons based upon three independent geometric pa- 
rameters (ki, K2, (f>). Continuously varying these parameters allows for a smooth 
transition between the different shape classes. This model employs continuum elas- 
ticity, differential geometry, and stationarity principles to predict how a naturally 
flat ribbon will deform when subject to arbitrary surface stress and/or internal 
residual strains. It establishes a quantitative relationship between surface stress, 
residual strain, elasticity, and helical properties such as helix angle, radius, axis, 
and chirality. The morphological transition from cylindrical helical ribbons to purely 
twisted ribbons and tubular structures can be tuned by relevant geometric param- 
eters. The predictions were validated experimentally with composites formed by 
bonding a pre-stretched sheet (or two sheets) of elastomer to an elastic layer of soft 
acrylic. The current results establish a new theoretical framework for predicting 
and tuning the shape of stress-driven helical structures for use in a broad spectrum 
of biological and engineering applications. 

6. Appendix: Mathematical description of the ribbon surface 

A mapping of an initially straight ribbon into the as-deformed, helical one, can 
be constructed as follows. First, recall that the centerline coordinates are given in 
Eq. (|2.14p . Now, consider a material point along the short edge of the ribbon in 
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the undeformed configuration given by Q°(s = 0,*) = P(0)+tB(0). Repeating the 
derivation in Eqs. (|2.3D through (|2.14l) in terms of t, where —w/2 < t < w/2, yields 
the following expression for the edge space curve upon deformation: 

Q(s = 0, t) = X (t)E x + Y (t)E y + Z (t)E z , (6.1) 

where 

Xo(t) = —z(od — sinat)(/ti — K2) sin cos 

B 2 

Y (t) = t -(erf — sin erf) 

Z (t) = 4 (cos at -1), (6.2) 

with a and B given in Eqs. (|2.9I) and (|2.10l) . respectively. Now, the ribbon surface 
can be constructed by "attaching" the space curve described by Eqs. (|A 1[) and (|A 2\ 
to each point s along the ribbon centerline, and locally rotating the space curve 
as dictated by the Frenet-Serret frame {T(s), N(s), B(s)}. That is, material points 
Q,°(s,t) in the undeformed configuration are mapped to Q(s, t) = X(s, t)E x + 
Y(s,t)~E y + Z(s,t)~E z upon deformation, such that 

Q(s,t) = P(s) + [X Q (t)E x +Y (t)E y + Z (t)-E z }A(s) 

= X{s,t)E x + Y(s,t)E y +Z{s,t)E z , (6.3) 
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where the rotation matrix A(s) is given by 

A(s) =T x (s)-E x ® E x + T y (s)E x <g> E y + T z (s)E x ® E z 

+ S a (s)E y ® E x + B y (s)E y ®E y + B z (s)E y <g> E z 
+ iV a (s)E 2 ® E x + N y (s)E z ®E y + N z (s)E z g> E 2 



/3 2 

1 ~(1 — cos as) 

a z 



E x ® E x + — ^(1 - cosas)E;E ® E„ - — sinasE^ <g> E z 
a^ a 



Bt 

+ -^(1 - cosas)E v <g> E s 



/? 2 



(1 — cos as) + cos as 



E y (g> E y H — sin asE y ® E z 



/3 T 

— sin asE z eg) E^ sin asE z ® E„ + cos asE z (g> E z 



(6.4) 



Finally, carrying out the matrix multiplications yields the following explicit expres- 
sions for X(s,t), Y(s,t) and Z(s,t): 



Bt 

X(s,t) = Xq(s) H t (at — sin at) 

a 13 



2/3 2 

1 H ^-(cosas — 1) 

a z 



( 1 — cos as)rt (1 — cos at) sin as 

a 

/3 2 

F(s,i) = ^b(s) H r[(r 2 - /3 2 )(1 - cosas) - a 2 cosas](at - sinai) 



/3 2 



(1 — cos as) + cos as 

8 



— ( 1 — cos at) sin as 



Z(s,t) = Zq(s) H (cos at — 1) cosas 



-[— 28 2 (at — sin at) + a 3 t] sin as . 



(6.5) 



where X (s), Y (s), and Z (s) are given in Eq. (|2.14D . 

Appendix A. Mathematical description of the ribbon surface 

A mapping of an initially straight ribbon into the as-deformed, helical one, can 
be constructed as follows. First, recall that the centerline coordinates are given in 
Eq. (|2.14p . Now, consider a material point along the short edge of the ribbon in 
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the undeformed configuration given by Q°(s = 0,*) = P(0)+tB(0). Repeating the 
derivation in Eqs. (|2.3D through (|2.14l) in terms of t, where —w/2 < t < w/2, yields 
the following expression for the edge space curve upon deformation: 

Q(s = 0, t) = X (t)E x + Y (t)E y + Z (t)E z , (A 1) 

where 

Xo(t) = —z(od — sinat)(/ti — K2) sin cos 

B 2 

Y (t) = t -(erf — sin erf) 

or 

Z (t) = 4 (cos at -1), (A 2) 

with a and B given in Eqs. (|2.9I) and (|2.10l) . respectively. Now, the ribbon surface 
can be constructed by "attaching" the space curve described by Eqs. (|A 1[) and (|A 2\ 
to each point s along the ribbon centerline, and locally rotating the space curve 
as dictated by the Frenet-Serret frame {T(s), N(s), B(s)}. That is, material points 
Q,°(s,t) in the undeformed configuration are mapped to Q(s, t) = X(s, t)E x + 
Y(s,t)~E y + Z(s,t)~E z upon deformation, such that 

Q(s,t) = P(s) + [X Q (t)E x +Y (t)E y + Z (t)-E z }A(s) 

= X{s,t)E x + Y(s,t)E y +Z{s,t)E z , (A3) 
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where the rotation matrix A(s) is given by 

A(s) =T x {s)-E x <g> E x + T y (s)~E x ®E y + T z (s)E x <g> E z 

+ B y (s)E y <g> E x + B y (s)Ej, <g) E y + B z (s)E a <g> E z 
+ iV K (s)E 2 ® E x + AT y (s)E, ® E y + A z (s)E z <g> E 2 



/3 2 

1 -(I — cos as) 

a z 

Bt 

+ ^(1 - cosas)E y (g> E x + 



E x (g) E x + — ?j-fl — cosaslE^ ® E„ — — sinasE^ ® E 2 
a 2 a 



(1 — cos as) + cos as 



Ey ® Ej, 



/3 T 

— sin asE z €3 E x sin asE z ® E„ + cos asE z (g) E z 



■ sinasE,, (gi E z 



(A 4) 



Finally, carrying out the matrix multiplications yields the following explicit expres- 
sions for X(s,t), Y(s,t) and Z(s,t): 



Bt 

X(s,t) — Xq[s) H -(at — sinai) 



2/3 2 

1 H ^(cosas — 1) 

or 



/3, 



(1 — cos as)rt (1 — cos at) sin as 

a 



B 2 

Y{s,t) = Y (s)-\ -Ut 2 - B 2 )(l - cosas) - a 2 cosas](at - sinat) 



+ t 



— x (1 — cos as) + cos as 
a z 





— r(l — cos at) sin as 

a J 



Z(s,t) = Z (s) H — - (cos at — 1) cosas 



-[— 2/3 2 (at — sin at) + a 3 t] sin as . 



(A 5) 



where A (s), io(s), and Z (s) are given in Eq. (|2.14l) . 
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